Hydrodynamic induced deformation and orientation of a microscopic elastic filament 



in 
o 
o 

(N 



M. Cosentino Lagomarsino, 1 I. Pagonabarraga, 2 and CP. Lowe 3 

'UMR 168 / Institut Curie, 26 rue d'Ulm 75005 Paris, Franc^ 
2 Departament de Fi'sica Fonamental, Universitat de Barcelona, C. Marti i Franques 1, 08028 Barcelona, Spaii^\ 
3 van 't Hoff Institute for Molecular Science, University of Amsterdam, 
Nieuwe Achtergracht 188, 1018 WV Amsterdam, The Netherlands 
(Dated: February 2, 2008) 

We describe simulations of a microscopic elastic filament immersed in a fluid and subject to a uniform external 
force. Our method accounts for the hydrodynamic coupling between the flow generated by the filament and the 
friction force it experiences. While models that neglect this coupling predict a drift in a straight configuration, 
our findings are very different. Notably, a force with a component perpendicular to the filament axis induces 
bending and perpendicular alignment. Moreover, with increasing force we observe four shape regimes, ranging 
from slight distortion to a state of tumbling motion that lacks a steady state. We also identify the appearance of 
marginally stable structures. Both the instability of these shapes and the observed alignment can be explained by 
the combined action of induced bending and non-local hydrodynamic interactions. Most of these effects should 
be experimentally relevant for stiff micro-filaments, such as microtubules. 
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PACS numbers: 87.16.Ka,47.15.Gf,46.32.+x,05.45.-a 

Semi-flexible polymers and filaments are important compo- 
nents of biological systems. All cytoskeletal filaments, used 
by cells for transport, morphology and force generation, fall 
into this category 1 1]. So do their assemblies, cilia and flag- 
ella, used to generate propulsion [2, 3]. Recent developments 
in experimental techniques allow the controlled synthesis, ma- 
nipulation and direct visualization of both real 0, S 01 and 
model filaments |7] of this type. Consequently, there is re- 
newed interest in the theory and simulation of both their static 
and dynamic properties, the structure and dynamics of DNA in 
flows [5] for example. If one is interested in dynamics, as we 
are here, one must consider that these filaments are normally 
suspended in a fluid. On the micron scale, inertia is generally 
negligible; any motion takes place in an environment effec- 
tively a billion times stickier than we experience in our daily 
lives . The dynamic response is determined by both elastic 
and fluid forces. Analytically, this is a non-trivial nonlinear 
problem. Nonetheless, at the linearized level the case of an 
elastic filament waved at one end has been solved |j4], and the 
predictions compared successfully with micro-manipulation 
experiments on an actin filament. The fluid was accounted for 
by introducing friction coefficients for parallel and perpendic- 
ular filament motion, independent of both the location along 
the filament and its configuration. Such an approximation is 
termed resistive force theory ||2|] and is satisfactory for several 
problems. Given experimentally observed flagella waveforms 
it gives good predictions for the swimming speed of sperma- 
tozoa 1 3 ] . It also predicts the buckling instability 1 9i] observed 
in sheared suspensions of filaments. 

However, resistive force theory is an approximation. In re- 
ality, a moving filament sets up spatially varying flow fields 
in the fluid that couple back to the motion of the filament it- 
self. To see why this coupling need not be trivial, suppose 
we approximate a filament of length L as a set of n beads 
separated by a fixed distance b = L/(n — 1), in the spirit of 
the "shish kebab" model [10] of a cylinder. A bead moving 



with velocity v experiences a hydrodynamic frictional force 
Fh = — 7o( v — v h), where 70 is a bead friction coefficient 
and vh is the velocity of the fluid at the location of the bead. 
Writing vh in terms of the flow fields of the other beads gives 
the hydrodynamic force on bead i 
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where is the vector connecting beads i and j 
hydrodynamic force acting on bead j and 77 the viscosity of 
the fluid Hill . Consider now a filament subject to a uniform 
external force density F x directed perpendicular to the axis of 
the rod (all forces and velocities are in the plane defined by 
the rod and the external force). If b << L equation Q gives 
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where s G [0, L] is the distance along the filament. The sec- 
ond term in the bracket is due to the non-local flow created by 
the rest of the beads. For a rigid body, the velocity is a con- 
stant, while the hydrodynamic friction force is greater towards 
the ends of the rod than in the middle. Thus, one expects a 
flexible filament to bend. In this letter, we examine this prob- 
lem numerically and show that this spatially varying interac- 
tion between the filament and the fluid substantially enriches 
the dynamic behavior of the system. Notably, the simulations 
predict hydrodynamic induced distortion and orientation that 
should be observable experimentally. 

We describe results from a simple numerical model that 
takes into account the interaction of the filament with its sur- 
rounding fluid more realistically than the resistive force ap- 
proximation. The position-dependent frictional forces are cal- 
culated explicitly from the flow fields generated by the fila- 
ment and will generally be configuration dependent and non- 
uniform. The semiflexible filament is treated as a curve, with 
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FIG. 1: Bending amplitude A, as a function of the the dimensionless 
force B, for a filament with L/b = 30 (solid line) and L/b = 100 
(dashed line). The insets are the corresponding shapes (the external 
force acts in the downward direction). 

its instantaneous configuration specified by a position vector 
r(s). The length is fixed and the elastic energy character- 
ized by a Hamiltonian H = n f Q C(s) 2 ds/2 111 211 . where 
C(s) = \d 2 r/ds 2 \ is the local curvature and k the stiff- 
ness. Numerically, the filament is modeled as a set of rigidly 
connected beads. The force acting on a bead is the sum of 
the external, tension, bending and hydrodynamic force. We 
calculate the latter from eq. Q, with Fj being the sum of 
all non-hydrodynamic forces acting on bead j. We can use 
eq. (|2j to calculate the friction coefficients of a rigid fila- 
ment modeled to this degree of approximation, yielding the 
friction coefficient 7 for motion perpendicular to the rod 
j x = 47rryi/ln (L/b), in the limit b « L [13]. If the bead 
spacing is interpreted as the cylinder radius this result is the 
same as the exact slender body hydrodynamic theory for a 
cylinder 0- An advantage of our approach is that, despite 
its simplicity, with increasing number of beads it approaches 
the correct result for a slender body, including the non-local 
nature of the hydrodynamic forces. It is accurate for very slen- 
der filaments and most of the biofilaments we are interested 
in satisfy this condition. Although a mathematically more 
sophisticated approach to treat the slender limit exists flUl . 
it does not account for the fact that the hydrodynamic force 
for a filament under a uniform field need not be uniform. As 
the force on a particle due its hydrodynamic interaction with 
its neighbors depends only on the external force and the in- 
stantaneous configuration of the filament, the simple implicit 
method described in ref. 1 3] suffices to integrate the equations 
of motion. These are solved subject to the condition that fila- 
ment inertial effects are irrelevant (the motion is overdamped, 
consistent with neglecting the fluid inertia). 

We first consider the motion of a filament under the action 
of a uniform field F K , acting perpendicular to the axis of an 
initially straight filament. The filament's shape evolves in time 
until a steady state is reached, where it drifts at a constant ve- 
locity with a fixed shape. That is, the bending, external and 
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FIG. 2: The perpendicular friction coefficient, ■y ± , relative to its 
value in the stiff limit 70 , as a function of the dimensionless force 
B. 



tension forces are balanced by the configuration dependent 
fluid force at every point along the filament. This steady state 
is a function of a dimensionless force B = L 3 F x /k. When 
B « 1, we expect the filament to behave as a rigid rod. 
With increasing B, significant bending will be required to bal- 
ance any non-uniformity in the hydrodynamic force. Typical 
shapes we observe for the steady state over a range of values 
of B are shown inset in Fig.^ In the steady state the filament 
is bent, indicating that the higher frictional force towards the 
end of the filament has to be balanced by a bending force. In 
the figure we plot a characteristic transverse distortion, A, de- 
fined as the distance between the uppermost and lowermost 
point of the filament along the direction of the applied force. 
We can identify four distinct regimes. For forces correspond- 
ing to small values of B (B < 50), the degree of distortion 
is small and increasing in proportion to B. This is the linear 
regime, where the bending can be understood as one of the 
lowest elastic modes H. II6II excited in response to the applied 
force (coming both from the external field and hydrodynamic 
interactions). For B > 100 the coupling between hydrody- 
namic and elastic forces is clearly non-linear, as evidenced by 
the rounded "U" shape of the filament in this regime. The 
bending amplitude saturates and, with increasing B, the fila- 
ment "U" shape becomes increasingly rounded. In this regime 
we also see that the friction coefficient starts to markedly de- 
crease, as shown in Fig. |2] This is because the progressive 
alignment of the filament with the applied force leads to an 
increasing fraction of parallel motion (see Fig.Q, character- 
ized by a lower friction. The perpendicular friction coeffi- 
cient is always greater than the parallel friction coefficient (not 
shown) but, for large B, appears to approach it I17II . Thus, in 
this limit the filament approaches hydrodynamically isotropic 
behavior. For even higher values of B, (B > 2000), we see 
yet different behavior. The filament initially adopts a new "W" 
shape, that, incidentally, has a more uniform hydrodynamic 
force density (and hence is closer to resistive force theory) 
and a lower distortion (shown by the lower data points in the 
figure) than its final state. The W configuration is marginally 
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stable; for any initial perturbation, after a transient time the W 
configuration rotates with the same handedness as the original 
perturbation and finally adopts a highly distorted but stable 
"horse shoe" shape (the upper data points shown in Fig. [Q. 
At still higher B (> 4000), we observe a regime (not shown) 
where following the formation of the "W" the filament again 
rotates but never adopts a new stable state. Rather, it exhibits 
a periodic zig-zagging motion indicating that above some crit- 
ical value of B there is no dynamically stable steady state that 
the filament can reach 1 18]. 

Now we consider a filament that is initially straight but with 
its axis tilted with respect to the force. If the hydrodynamic 
friction along the rod is uniform, the filament will remain 
straight and maintain the same orientation. It will move at 
the constant velocity for which the friction and external forces 
balance each other 1 19]. Because the parallel and perpendicu- 
lar friction coefficient differ, this will generally be at an angle 
to the external force. A truly rigid rod, even with non-local 
hydrodynamic interactions, must also move at constant speed 
maintaining its initial orientation, as is the case for any rigid 
object. However, as shown in Fig. [3] we find that if the ini- 
tial configuration is rotated anti-clockwise with respect to the 
external force, as the rod translates it rotates clockwise until 
its ends are again aligned perpendicular to the force. During 
the rotation the filament exhibits a transient drift motion in a 
direction inclined to the force but in the steady state, unlike 
a rigid rod, the drift is along the applied force direction. We 
recover identically the case discussed above. 

The orientation angle as a function of time - after a time dur- 
ing which bending is established - decays exponentially 1 20 ] . 
The characteristic time for this decay, th, relative to the time a 
rigid rod takes to translate its length (tt = 7 ± /-F x ) is plotted 
as a function of B in Fig. [3] 

This effect is related to the flexibility, because a truly rigid 
rod cannot display this behavior. Why should a flexible fila- 
ment behave so differently? As long as the rod is not aligned 
parallel to the applied force, there is a component of the force 
perpendicular to the filament inducing bending. Taking the 
situation illustrated in Fig. [3] the bending will slightly align 
the tangent to the filament at the right hand end parallel to the 
applied force and the tangent at the left hand end perpendic- 
ular. The local friction coefficient is higher perpendicular to 
the filament than parallel. So, even if the two ends move with 
the same velocity, the drag force on the right hand end will 
be slightly lower that that on the left. Thus, a torque tends 
to rotate the filament clockwise, as we observe (see Fig. QJ. 
A simple semi-quantitative analysis indicates that in the lin- 
ear regime, where the degree of bending is proportional to B, 
the torque generated is proportional to B. One therefore ex- 
pects that the time it takes the filament to reorient scales as 
1/B, as shown in Fig. [3] Note that this simple argument only 
invokes the non-local hydrodynamics to generate the bend- 
ing that breaks the symmetry between the two ends. Arguing 
along the same lines, resistive force level hydrodynamics pre- 
dicts rotation of a bent filament. It fails because it does not 
include the mechanism that generates the bending. 
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FIG. 3: Log-log plot of the characteristic hydrodynamic re- 
orientation time relative to the translational time, th, as a function of 
the dimensionless force B. Inset are the corresponding motions for 
two characteristic values of B, viewed in the center of mass frame of 
reference, the initial and final states being I and F respectively (the 
external force acts in the downward direction). 



This analysis has some interesting consequences. It appears 
that any degree of elasticity will induce the rotation of the fil- 
ament. Hence, the completely rigid limit is singular, in the 
sense that the absence of orientation is because the reorienta- 
tion time becomes indefinitely long as the rigidity increases. 
It also implies that the motion of a flexible filament aligned 
parallel to the force will be unstable. A small perturbation of 
the angle from the parallel, with a certain handedness, gener- 
ates a torque with the same handedness. The perturbation will 
be amplified and the filament will rotate until the dynamically 
stable perpendicularly aligned "U" state is reached. This rea- 
soning also explains why the "W" shape is marginally stable. 
By symmetry, the two outer "U" sections cannot contribute 
any net torque. However, the central convex section will dis- 
play an instability with respect to any perturbation from the 
perpendicular, consistent with our observation that the "W" 
configuration eventually rotates to form a horse shoe. 

Could this behavior be relevant in practice? Appreciable 
bending requires B » 1 while hydrodynamic induced ro- 
tation should occur for any value of B (although for small B 
the reorientation time will be long). In the case of biofila- 
ments, a gravitational field only gives low values of B. Mi- 
crotubules have a stiffness of k ~ 50pN/jm 2 l2lll . so for 
a typical length (10/im) we estimate that B < 10~ 3 . Ac- 
celerations three orders of magnitude greater than gravity are 
needed to reach even B ~ 1 (which can be achieved in a cen- 
trifuge). However, bio-filaments are commonly charged, and 
hence will react to applied electric fields. If we concentrate 
on microtubules, they are negatively charged with an effec- 
tive charge density q of approximately 200e//im |6], where 
e is the elementary charge. The strongest electric fields for 
which microtubules are stable, reported in ref [6], are of order 
10 4 V/ m, corresponding to a force density of 0.3pN/ [i,m. For 
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a microtubule of length L ~ h^m (as used in the experiments 
reported in 1 6]) one then finds B ~ 1. However, since B oc L 3 
for longer microtubules of L ~ 30/jto, we have _B ~ 200. It is 
therefore possible to apply forces that should induce observ- 
able bending. A further condition to observe the orientational 
behavior we describe is that the hydrodynamic orientational 
time th is shorter than any other relaxation time in the prob- 
lem. For example, we have thus far ignored rotational diffu- 
sion. It will tend to randomize the orientation on a charac- 
teristic time scale te> ~ 7 L 2 /kT. From Fig.[3]we see that 
th ~ 7 ± /(F x B), so t h /t d ~ L/{B 2 X), where A(= n/kT) 
is the persistence length. One can therefore neglect diffusion 
if L/(B 2 \) « 1. For a stiff filament, L/X « 1 (for a 
lOfim microtubule this ratio is 10~ 2 ), so this condition is au- 
tomatically satisfied if B > 1. In addition, microtubules have 
an electric dipole. This will lead to a torque tending to align 
them parallel to the electric field while they translate and bend. 
The ratio of the dipole re-orientation time Td to the transla- 
tional time is t^/tt = qL 2 /d, where d is the electric dipole. 
Experimental results indicate that d ~ eL 0]. Hence the ra- 
tio is Td/rx ~ 200L/Ad with the length Ad = l/.im. So for a 
microtubule with L > lfim the dipolar re-orientation should 
be negligible compared to hydrodynamic orientation. 

An analysis along these lines suggests that for actin and 
DNA it is possible to achieve B » 1. However, the con- 
dition that the hydrodynamic orientation time is shorter than 
all other time-scales is not easily satisfied. There will be a 
competition between the effects we describe here, acting to 
distort and align the filament, and thermal effects, acting to 
randomize the orientation and maintain the equilibrium struc- 
ture. Finally, the situation with carbon nanotubes is more flex- 
ible, given the greater control of physical properties of these 
objects 1 22]. It may be possible to reach the B >> 1 regime 
where the instability sets in. 

In this letter we have shown how the non-local nature of the 
hydrodynamic interactions affects the dynamics of inextensi- 
ble elastic filaments subject to a uniform external field. The 
method captures all the relevant configurational couplings. 
Although the hydrodynamic treatment is only exact for an in- 
finitely slender rod, it is computationally simple compared to 
other techniques 1151 12311 . The non-uniform, configuration- 
dependent friction the fluid exerts on the filament induces dis- 
tortion and a corresponding increase in mobility. In the dy- 
namic steady state, the degree of bending depends on the stiff- 
ness of the filament, a fact that could be used experimentally 
to determine k. There is a crossover from the linear regime 
to a plateau value for the degree of distortion, reminiscent of 
other pattern-forming systems. Our numerical results suggest 
that at still higher degrees of forcing there exist long-lived 
marginally stable states and that eventually the filament be- 
haves in an unsteady but hydrodynamically isotropic manner. 
The fact that fluid friction induces bending means that a flex- 
ible rod aligns perpendicular to an applied force. This is not 
expected if one neglects either the flexibility of the filament or 
the non-local hydrodynamics. It does, however, offer a novel 
route by which one may manipulate the orientations of fila- 



ments experimentally. 
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